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We study diffusion-limited (on-site) pair annihilation A + 
A — > and (on-site) fusion A -\- A ^ A which we show 
to be equivalent for arbitrary space- dependent diffusion and 
reaction rates. For one-dimensional lattices with nearest 
neighbour hopping we find that in the limit of infinite re- 
action rate the time-dependent n-point density correlations 
for many-particle initial states are determined by the correla- 
tion functions of a dual diffusion-limited annihilation process 
with at most 2n particles initially. Furthermore, by refor- 
mulating general properties of annihilating random walks in 
one dimension in terms of fermionic anticommutation rela- 
tions we derive an exact representation for these correlation 
functions in terms of conditional probabilities for a single 
particle performing a random walk with dual hopping rates. 
This allows for the exact and explicit calculation of a wide 
range of universal and non-universal types of behaviour for 
the decay of the density and density correlations. 
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1 Introduction 



In this paper we investigate a class of models describing diffusion-limited 
annihilation of identical particles. This is a process where particles of a single 
species diffuse (rather than move ballistically) and undergo an annihilation 
reaction when two particles meet. The outcome of such a reaction may be 
either a single particle, or both particles disappear. This process can describe 
both chemical reactions where the particles change their state into an inert 
reaction product which takes no part in the subsequent dynamics of the 
system, or physical reactions where the particles actually annihilate under 
the emission of radiation. On an abstract level, both kinds of reactions are 
identical. A known one-dimensional process with "physical" annihilation is 
e.g. laser- induced exciton dynamics on polymers Clearly, one is not 
interested in the equilibrium behaviour of such a model, since in the absence 
of continued particle production all particles eventually disappear. Instead 
one would like to understand non-equilibrium properties such as the time- 
dependence of the particle concentration for a given initial condition and the 
degree of universality obtained from specific choices of models. 

Diffusion-limited pair annihilation of particles of a single species is a well- 
studied system (for a recent review see H), but there are still interesting 
open questions particularly in the presence of spatially inhomogeneous hop- 
ping and reaction rates which have so far not been addressed. In homo- 
geneous, translationally invariant environments the particle density decays 
with a power law which depends on the dimensionality of the system. Both 
theoretically and experimentally one finds p{t) ~ 1/^/Dt in one dimension. 
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Interestingly, this result is at variance with the dimensionality- mdependent 
mean-field behaviour p(t) ~ l/{Dt) which is correct only in three (and 
higher) dimensions. The amplitude of this decay is universal in the sense 
that it depends neither on the initial density for random initial conditions 
nor on the reaction rate 0. However, if particles are moving in an arbi- 
trary, non-translationally invariant energy landscape it is not obvious how 
this will change the decay of the local or overall particle concentration. In 
our study which is a continuation of previous work we first discuss uni- 
versality of a general class of models, but then shall pay special attention 
to one- dimensional systems with nearest neighbour interaction. The physi- 
cal motivation behind the study of one-dimensional system is not only their 
experimental relevance for polymer physics, but also their theoretical im- 
portance in the understanding of the role of fiuctuations in low-dimensional 
systems. In both one- and two-dimensional systems diffusive mixing is in- 
efficient and leads to the building up of large-scale correlations. Thus the 
classical mean-field rate equations for the study of these systems tend to fail 
and require a more sophisticated treatment. 

To this end we investigate by exact methods a family of lattice models 
with space-dependent hopping and pair annihilation rates, as one would have 
e.g. in disordered systems or in the presence of a space-dependent external 
forces affecting the diffusive motion of the particles. These particles have no 
attractive or repulsive interaction between themselves, they hop with fixed 
rates d^y from lattice site x to site y. When two particles meet on site x 
they both annihilate with a rate A^;. We address two different problems: (1) 
The first is to identify other systems which are in the same universality class 
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in any dimension. In fact, our approach of using similarity transformations 
extends to systems defined on arbritrary lattices. This idea is not new for 
systems with exclusion dynamics @] - [@|. Here we extend the application of 
similarity transformations to models without site-exclusion. This generaliza- 
tion is useful for a renormalization group treatment of the problem ||^, H, 10 



as it yields exact relations for expectation values for one model in terms of 
expectation values for the related system without any further RG analysis 
for each of these different models. (2) Since the RG treatment shows that 
the fixed point of the homogeneous system with nearest neighbour hopping is 
the limit of infinite reaction rate, and since mean-field is particularly poor in 
one dimension, the second problem we address is a more detailed investiga- 
tion of the one- dimensional case in the limit of infinite reaction rate but with 
(arbitrary) space-dependent nearest neighbour hopping rates. Here we use 
not only similarity transformations but also free fermion techniques which, 
as shown below, are a convenient framework to formulate general properties 
of annihilating random walks in one dimension. From our treatment we find 
that the time-dependent density correlations for an arbitrary many-particle 
initial state have an exact representation in terms of conditional probabilities 
for a single particle performing a random walk with dual hopping rates. In 
the case of constant hopping rates the predictions of this model are in ex- 
cellent agreement with experimental data on exciton dynamics on very long 
ordered polymer chains. Thus we expect that our model gives an equally 
good description of the behaviour of inhomogeneous systems. 

The paper is organized as follows. In Sec. 2 we briefly review the de- 
scription of reaction-diffusion processes in terms of a master equation which 



is written in a quantum Hamiltonian formalism. In Sec. 3 equivalences be- 
tween models without site exclusion are studied. This is extended in Sec. 

4 to the limit of infinite reaction rate and nearest neighbour hopping and 
the main results of the paper are derived. In Sec. 5 some specific results 
for random initial conditions are derived and in Sec. 6 the main results are 
summarized and discussed. 

2 Quantum Hamiltonian formalism for stochas- 
tic reaction-diffusion systems 

The reaction-diffusion processes considered in this work are systems of clas- 
sical interacting particles which hop stochastically on a lattice and undergo 
some reaction when they meet on the same lattice site. For definiteness we 
shall consider only systems of one species of particles even though most of 
what is discussed in this and the following Section can be generalized to sys- 
tems of many species. In symbolic presentation particles are denoted by A, 
a vacant site by 0. The expression kA denotes the presence of k particles are 
in the same site. 

We define the process in terms of a master equation for the probability 
P{rj;t) of finding, at time t, any configuration r) of particles on a lattice 

5 oi C sites. Here rj — {77(1), 77(2), 77(£)} where ri{x) are the integer- 
valued particle occupation numbers at site x. The lattice is at this stage 
arbitrary. Later, when discussing one-dimensional systems we consider only 
linear chains with C = L sites and periodic boundary conditions. For clarity 
we shall use in this case site labels k, I, m instead of x, y, z. 
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A general master equation reads 
d 



^^P{V'- t) = Y, {w{r]', v)P{v, t) - wiv, v')P{v'- 1)} (1) 



The positive contribution to the (infinitesimal) change in probability is the 
sum of probabilities that the system has been in a state i] times the rate 
w{ri',ri) of fiipping to r]'. The rates satisfy < w{rj',vi) < oo and define 
the exponentially distributed life time T~^{rj') = J2rif^ri' '^(VjV') of ^ state 
1]'. This is the total rate for the state rj' to fiip into any other configuration 
and hence represents the loss in probability for the state rj' in Eq. (|l]). We 
shall now express the time evolution given by the master equation in terms 
of a quantum Hamiltonian H. This is discussed in detail in and also in 



earlier work [0, |T5|, |T6[ and we shall repeat only the essential elements of the 
mapping. The advantage of this approach is that there are standard methods 
of dealing with the resulting time evolution operator H. The applicability 
of these techniques, in the case at hand primarily similarity transformations 
and the renormalization group technique for systems without site exclusion 
and free fermion techniques for models with site exclusion, does not arise 
naturally if the master equation is written down in the standard form ([^). 

2.1 Definitions 

The idea is to represent each of the possible particle configurations t] in the 
set X = by a vector | rj ) which together with the transposed vectors 
(r/l form an orthonormal basis of a vector space X = (C°°)®^. The basis 
is chosen such that for a single site the state |n) with ri{x) = n particles 
is represented by the unit column vector with a 1 at the n*^ position and 



zero elsewhere. Then the state of the lattice as a whole is represented by 
the tensor product Itj) = \ r]{l) ) . . . ® | f]{C) ). Now one can conveniently 
represent the probability distribution by a state vector 

\P{t)) = J2P(v;t)\v)- (2) 

riex 

Using P{ri; t) = {ri\ P(t) ) one writes the master equation (|I]) in the form 

^^P{v;t) = -{v\H\P{t)) (3) 

where the off-diagonal matrix elements of H are the (negative) transition 
rates and the diagonal entries are the life times of the states. 

Therefore the time evolution of this probability vector is given in terms 
of a linear 'Hamilton' operator H acting on X 

jy{t)) = -H\P{t)). (4) 

A state at time t = + t is given in terms of an initial state at time to by 

|P(to + r))=e-^-|P(to)). (5) 
Note that ( s \ P{t) ) = ^^^^ P(r/; t) = l\/t where 



s 



vex 



This relation expresses conservation of probability and implies 

{s\H = (7) 

for any stochastic process. We shall call a matrix which satisfies (|^) and 
Hn,r]' < for all off-diagonal elements a stochastic Hamiltonian. 
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This definition is manifestly basis dependent. If a stocliastic Hamiltonian 
H is related to some other Hamiltonian H by a, similarity transformation, 

H = BHB'^ (8) 

we shall call the systems equivalent, irrespective of whether H is stochastic 
or not. If H is stochastic, then the two processes are called equivalent. An 
equally useful relation is 

H = BH^B-^ (9) 

relating a stochastic process to the transpose of a matrix H. If both H and 
H are stochastic, we shall say that these two processes are enantiodromic 
with respect to each other. Any family of transformations B generates a 
family of related processes (equivalent or enantiodromic). Any one member 
of this family is called a representative. 

Expectation values {Q{t)) = {s \Q\ P{t)) are calculated as matrix ele- 
ments of suitably chosen diagonal operators Q. A complete set of observables 
are the occupation numbers r){x). Defining projection operators on states 
with n particle on site x of the lattice as 

QM-\n)AnL (10) 

one finds that the density Px{t) of particles at site x is given by the expec- 
tation value (rix) of the operator rix — Yl'^=o''^Qxi'>T)- This is the diagonal 
matrix n acting non-trivially only on site x with matrix elements {n)ki = kSk,i- 
Density orrelation functions ( ■ ■ ■ Ux^ ) are computed analogously. 

For later convenience we also introduce the operators a^. The single-site 
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matrix a"*" has elements {a'^)m,n = ^m.n+i- Hence 

a^|n) = |n+l), {n\a~^ = {n — 1\. (11) 

creates a particle with unit rate when acting to the right. The single-site 
matrix a~ has elements {a~)m,n = ^Sm,n-i- Thus 

a^\n) = n\n — 1) , {n\a'~ = {n + l){n + 1\. (12) 

annihilates a particle, but with a rate proportional to the occupation number 
prior to the annihilation. This corresponds to annihilating any one of the 
existing identical particles with unit rate. It is understood in these equations 
that a~\0) = (Ola"*" = 0. Note that the matrices satisfy a harmonic 
oscillator algebra a~a+ — a+a~ = 1. The product n — a+a~ is the number 
operator. Matrices acting on different sites commute. 

The state ( s | is the coherent state {s\ = (0|exp(A~) where ~ 
Ylx^t ^^'^ (0| = ^he row vector representing the completely 

empty lattice. Using the harmonic oscillator algebra and the factorization of 
exp {A~) one finds the following relations 

{s\a^ = {s\, {s\a- = {s\n^ (13) 

which are useful for the construction of the Hamiltonian H and for the deriva- 
tion of the equations of motion for correlation functions. 

In the infinite reaction limit discussed below, the system we investigate 
reduces to a two-state model where each lattice site is either empty or occu- 
pied and the state space reduces to X = (C^)®"^. For this situation is more 
convenient to use Pauli matrices cr^ acting on site k for the description of 
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the elementary events of the stochastic time evolution. A complete set of 
observables are now the occupation numbers Uk = 0,1. Projectors on states 
with a particle on site k of the chain are then defined by the diagonal matrix 

nk = l{l- crl) . (14) 

One may equally well use spin language. In this interpretation projects 
on a spin down at site k and Vk = 1 — rik projects on vacancies or spin up 
respectively. 

Changes in the occupation numbers are effected by = {a^ ±ial)/2. In 
our convention creates a particle at site k when acting to the right, while 
annihilates a particle at site /c.Q Note that 

(s|Sj^ = (s|nfc and ( s = ( s |(1 - n^). (15) 

Introducing the ladder operators = ^^^^ one may write 

(s| = (0|e^^ (16) 

Using the commutation relations for the Pauli matrices then yields (p!5|). 

Some particular initial states and initial distributions used below: An ini- 
tial state with particles place on sites xi, . . . , xat is denoted by | xi, . . . , xat ). 
An uncorrelated initial state with a site-independent Poisson distribution is 
denoted by |p) where p stands collectively for all average site-densities px- 
The vector representing this inhomogeneous product measure is the factor- 
ized state \ p) = Ylx I Poi^) ) where | Po{x) ) = exp [pxiat )• For two- 

state models the same notation is used for an uncorrelated initial state with 

^This somewhat counterintuitive notation results from long standing convention which 
we do not wish to change. 
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bimodal distribution. In this case, the probabihty of finding site x occupied 
is px and the probabihty of finding it vacant is 1 — Px- The completely empty 
lattice is always represented by the vector | ) . 

2.2 Construction of the quantum Hamiltonian 

Now we are in a position to derive H. The process we consider as a repre- 
sentative is defined by following elementary rules: 

• Particles move on a lattice from sites x ^ y with a fixed rate p{x, y) 
which is not necessarily symmetric. 

• Any two particles on a site can annihilate with rate \{x). 

The idea behind these rules is the description of particles which have no 
physical interaction (repulsive or attractive), but only a chemical interaction 
which allows the formation of an inert product from two single particles 
which then has no infiuence on the subsequent dynamics. In this model no 
assumptions are made about the physical or chemical origin of the space 
dependence of these rates. 

The precise nature of the dynamics is defined by a master equation 
With matrices defined above the stochastic quantum Hamiltonian for the 
process has the form = H[) + Hj^ where H[) describes hopping Hj^ = 
J2x describes annihilation. It is given by 

H^ = - p{x,y){a-a^ -Tlx) -^\{x) {{a-f -nxirix-l)) . (17) 

x^yeS xGS 

The first sum runs over all pairs {x, y) of sites of the lattice 5* and describes 
hopping of non-interacting particles from x to y whereas the second sum 
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runs over all sites of the lattice and gives the annihilation events. According 
to the general rules the off-diagonal part represents the possible moves the 
system can perform in an infinitesimal time step. The diagonal part ensures 
that probability is conserved. Using ( ]IB| ) and the harmonic oscillator algebra 
it is easy to check that (0) holds for each elementary process and hence for 
H. 

Consider now the infinite annihilation limit A(x) = A(x)/Ao with Aq — > 0. 
In practical terms this means that the time scale Aq on which the annihilation 
takes place is much shorter than the time scales p{x, y) set by the hopping 
events. In this limit two things happen: Firstly, in the initial state any 
site occupancy by an even number of particles is immediately reduced to a 
vacancy and any site with an odd number of particles is left with just a single 
particle. Secondly, if by a hopping event a double occupancy is created, this 
results immediately in a pair of empty sites. Thus an attempted hopping 
event of the form A^Ay — » ^x'^Ay with rate y) is equivalent to the pair 
annihilation process A^Ay ^x^y with rate p{x,y). In short, the system 
becomes equivalent to an exclusion process with particle hopping and pair 
annihilation with annihilation rates that are equal to the sum of the hopping 
rates p{x,y) +p{y,x). The Hamiltonian for this process reads in terms of 
Pauli matrices 

H'' = - Ep(^' y^^'^'y + '^4 - ^^)- (18) 

This can be derived in a more formal way using the infinite rate formalism of 



1^ . The two systems I and II defined by ([TtD , ([TsI) are subject of this paper 



The infinite rate model (0) will be studied in some detail in one dimension 
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with nearest neighbour events {p{x, y) = for y ^ x ±1). 

Since the process conserves particle number modulo 2, it splits into two 
distinct subsectors corresponding to even and odd numbers of particles. For 
the two state process (0) the operator 

Q = tl^k = {-if (19) 

fc=i 

has eigenvalue ±1 in the even (odd) sector respectively. Hence = (1 ± 
Q)/2 projects on these sectors. Corresponding to this separation of the 
dynamics it is convenient to use 

( s ={s\P^ (20) 

instead of ( s | for the calculation of expectation values. 

3 Similarity transformations for DLPA 

We want to consider equivalences between model I and other systems of par- 
ticles hopping on a lattice and which, independently of the hopping process, 
undergo some reaction when they meet on the same lattice site. I.e. we 
are looking for a stochastic process of the form H = Ho + Hr such that 
H = BHB~^ and where Ho describes hopping on non-interacting particles 
and where Hr = J2x describes some on-site chemical interaction. At- 
tempting to identify all possible equivalent systems appears to be a hopeless 
enterprise which is the reason why we restrict ourselves to equivalent pro- 
cesses with only on-site interaction. Therefore we assume B to factorize, 
B = i?i . . . ® Be- This has a consequence for the choice of itself: 
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The transformed hopping Hamiltonian Hd should not contain two-site op- 
erators other than the hopping matrices Qj^ Qjy SiS such objects not describe 
non-interacting particles or strictly local interactions as demanded. This 
requirement is met by a matrix — exp {a^a' + Pxa'^ + ^x^x + Sx) with 
arbitrary parameters Oix, (3x,lx,Sx- Since (s \ has to remain invariant under 
the transformation, one has to choose jx = —Q^x and 5x = —Px- This leaves 
a family of transformations with 2C free parameters. 
Using the harmonic oscillator algebra one shows that 

Ba-B-' = e"a- + -(l-e") (21) 

a 

Ba+B-^ = e-"a+ + l-e-°^ (22) 

which are necessary for the calculation of H. Applying the transformation 
to Hj^ shows that positivity and reality of the transformed rates constrains 
the transformation parameters to be in the domain — In 2 < a^. < 0, /S^; = 0. 
Any transformation with (3x^0 leads to events with negative rate. Even 
though this can be compensated by inclusion of particle creation events in 
the original Hamiltonian, we do not consider this possibility in the present 
work.^ li (3 — the hopping part Hj^ is invariant under the transformation if 
one chooses the transformation parameters ax = a space- independent. Thus 
a transformation meeting all the requirements discussed above gives rise to 
the one-parameter family of processes 

^ - ^ p{x,y){a^a^ -Ux) 

^We just note for illustration that inclusion of single particle annihilation ^ — > and 

pair creation 2A with suitably chosen rates is equivalent to a process with pair 
annihilation, branching A 2A and creation A under a transformation with /3 ^ 0. 
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- J2 K^) [(1 - & + bat)^)' - - 1)] • (23) 

where particle either annihilate 2 A with rate (1 — 6) or undergo a fu- 
sion reaction 2A — * A with rate bX^ where b = 2(1 — exp(«)). This shows that 
pair annihilation and fusion are in the same universality class irrespective of 
the lattice on which the system is defined and also irrespective of the physical 
environment specifying the rates p{x,y) and X{x). The choice of lattice is 
completely arbitrary. The equivalence of the two processes was first discov- 
ered by Krebs et al. ^ for a model with site-exclusion and homogeneous 
nearest neighbour hopping in one dimension. Universality for finite reaction 
rate in any dimension was shown by Peliti using the RG approach. 

Since the mapping defined by the transformation B is exact it is worth- 
while to investigate some of the consequences for expectation values for the 
two processes. First we note that both ( s \ and B factorize and therefore 
{s\B = {s\. To proceed, we consider as initial state uncorrelated random 
initial conditions with site-dependent density px, i.e. we assume that at time 
zero particles have been distributed randomly and independently at each site 
with a Poisson distribution Px{n) = exp (— p^)p"/n!. This can be achieved 
by filling an originally empty lattice by the process H = — J2x Px/'T'i.O't ^ 1) 
up to a time t = t. Then the creation is switched of and annihilation and 
fusion take place. For the pure annihilation process the concentration Py{t) 
of particles at site y is given by Py(t) = { s In^exp {—HH) \ p). In order to 
relate the density 

p,(t) = (sK5-V^'*S|p) (24) 
for the mixed process to Py{t) we have to study the effect of the transforma- 
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tion on both | p ) and Uy. A short calculation shows that B\ Pq{x) ) remains a 
Poisson distribution, but with density px = Px exp (—a). On the other hand, 
(|13D together with (^TJ) yields {s\BnyB^^ = {s\a~B^^ = exp{a){s\ny. 
This gives 

Pyit) = e^yit) (25) 

where the initial densities for the mixed process are related to the initial 
density of the pure annihilation process by the relation p^ = p2^exp(— a). 
For a A;-point correlation function F{xi, . . . ,Xk',t) = { rix^ (t) ■ ■ ■ (^) ) where 
all Xi are pairwise different, (^) generalizes to 

=e'^"F(xi,...,Xfc;t). (26) 

If two or more coordinates are equal, one has to use (^Tj) and (^21). For the 
two-point function F2{y; t) = { riyit) ) this yields 

F2iy; t) = e"^F2iy; t) - e"(l - e")p,(t). (27) 

It is interesting to consider the ratio A{t) = {N^{t)) - {N{t)f/{N{t)) 
of the particle number fluctuations to the particle number ( ) = X] ( ) 
at time t. For the mixed process one obtains the exact relation A(t) = 
exp (a)(A(t) — 1) + 1. This relation can be used to measure the branching 
ratio of an experimental system, if one can measure both the average particle 
number and its fluctuations and if A for the pure annihilation process can be 
calculated analytically for the respective physical environment. In the sim- 
plest case of infinite pure annihilation and homogeneous nearest neighbour 
hopping in one dimension this quantity can be calculated exactly |TB| and 
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is for sufficiently large times independent of time and of the initial density, 
A(t) ^2-^2. 

4 Infinite reaction limit 

The transformation of the last section is independent of X{x). Therefore 
all the results translate with little modification into the dynamics of the 
exclusion process H^^ defined in (|T8|) which describes the dynamics with 
infinite reaction rate. One just has to replace in the transformation 

B, the results concerning the relations between correlation functions remain 
unchanged. 

However, this simple observation is not all there is to say about this 
limit. It is known that in one dimension with nearest neighbour hopping 
DLPA is related to zero-temperature Glauber dynamics by a domain-wall 



duality transformation ^7j. This transformation was shown recently to be an 
invertible similarity transformation |[T^ . On the other hand, it was noted that 
zero-temperature Glauber dynamics can be brought by another similarity 
transformation into a form which is the transpose of the Hamiltonian for 
DLPA 0. From this we conclude that there must be a matrix V such that 

H'' = V-\H'YV, (28) 

i.e. DLPA in one dimension with nearest neighbour hopping is self-enantiodromic 
in addition to being equivalent to the mixed pair- annihilation/fusion pro- 
cess.0 

^ Strictly speaking, this applies only to the sector with an even number of particles. 
For an odd number of particles the situation is slightly more complicated. 
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From now on we consider only the one- dimensional process on a ring with 
L sites and with arbitrary nearest neighbour hopping rates Ik = p{k, k — 1) 
and Tk = p{k, A; + 1). In order to avoid unnecessary technical complications 
with boundary terms, we consider only the sector with an even number of 
particles. Combining the results of Refs. [|l^, 0] yields 

V = 7i72 . . . 72L-1 (29) 

where 

72.-1 = ^[{l + t)al-il-t)] (30) 

72. = + (31) 

To prove the enantiodromy relation (|28|) one notes that V is unitary and 
transforms Pauli matrices as follows: 

= {^ij (32) 

where Q = Y[k=i ^1 — For the even particle sector where Q = 1 the 

relation ( pSD follows for constant rates Ik = Vk = D. 

In order to apply enantiodromy to inhomogeneous systems we first rewrite 
the Hamiltonian ( p!8[ ) adapted to the present case. For hopping rates as 
defined above the Hamiltonian reads 

L 

H = -Y.^rkhl + hK) (34) 

k=l 
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with the hopping-annihilation matrix /i^ = s~l^s^_^^ + s^s^_^-^ — Uk- Applying 
the transformation yields the process H = V^^H^V given by 

L 

H = -Y.{nK + khU). (35) 

k=l 

The enantiodromic process is of the same form ( P^ as the original process, 
but with dual hopping rates 

^k = rk , rk = h+i- (36) 

We shall refer to the environment defined by the dual rates as to the dual 
environment. The sector with an odd number of particles can be transformed 
in a similar way using V = T>a^. Some care needs to be taken in the 
treatment of the boundary term. 

In order to make practical use of the enantiodromy we note 

{sr-n,, . . .n,„e-^*| Po ) = ^{Po l^^e-^*(l-<-iO • • • {l-at^.,atjV-'\ s) 

(37) 

This may be written in the form 

( s r^n,, . . . n,^e-^'\ Po) = {s \Qe-^'\ ) (38) 

where Q is determined by the original initial condition and the initial state 
\Pq) is determined by the set of sites {ki, . . . , km}- The importance of this 
result is seen in the transformation law s )^^*^" = — ). Hence 

I Pg ) is a linear combination of initial states of at most 2m particles. Since 
H does not have any particle creation terms the time- dependence of the 
m-point correlation function with an arbitrary many-particle initial state is 
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completely determined by the dynamics of the system with not more than 
2m particles. In particular, the density at site k is given by the dynamics of 
the dual system with just two particles placed initially at neighbouring sites 
k-l,k. 

It remains to consider the dynamics of Q. Without specifying the original 
initial distribution | Pq ) and using (|T5|) we can always write Q as a linear 
combination of products with an even number of particle annihilation oper- 
ators . . . where 2p < 2m since H does not create particles. Now one 
can follow standard procedure and perform a Jordan- Wigner transformation 
^ by introducing Qk = Y[i=i '^i fermionic annihilation and 



creation operators 

4 = Qfc-i , Cfc = Qk-ist (39) 

satisfying the anticommutation relations {c^, q} = {c|,, c]} = and {cj,, q} = 
6k^i- In terms of these operators Q can be written as a linear combination of 
products with an even number of fermionic particle annihilation operators 
. . . and the initial state is of the form c|,^ . . . cl^^ \ ). The Hamiltonian 
H is bilinear in ct, Ck and leads as shown in Ml to a linear time-evolution 



equation for Ck{t) = exp {Ht)ckexp {—Ht). The solution of the differential- 
difference equation obtained by taking the time-derivative of Ck{t) is a single- 
particle problem, viz. the solution of the initial value problem of a single 
random walker in the hopping environment defined by H, i.e. 

Ckit) =J2Pik;t\l,0)ci (40) 

I 

where P{k] t\l, 0) = {k\ exp {—Ht) \ I ) is the conditional probability for the 
single-particle problem. This reduces the calculation of the expectation value 

19 



of Q to the calculation of correlators of the form ( s |c/^ . . . ci^pC\^ ■ ■ ■ c^.^^ | ) at 
time t — Q which are given by the anticommutation relations and c;| 0) = 0. 

The appearance of free fermions in this problem of stochastic dynamics of 
classical interacting particles may seem surprising. However, to calculate the 
two-particle matrix transition probability {m,n |e~^*| k, I ) one can either use 
the free fermion description, or, in a less technical way, remind oneselves of 
the meaning of an annihilating random walk and the description of random 
walks in terms of a sum over the canonical path space. In discrete space 
and time the transition probability (or conditional probability) for a single 
particle P{m;t\k,0) is the sum over all paths leading from k to m, each 
weighted with its proper statistical weight given by the hopping rates and the 
particular form of the trajectory. If two non-interacting particles, one starting 
at site k and the other at site I, move, then the transition probability that 
the particle which started at site k < I reaches site m < n and the particles 
which started at site / reaches site n at time t is still the sum over all possible 
trajectories which connect k with m and / with n, where each single trajectory 
has the same weight as in the single particle case. Hence, for non-interacting 
particles, P{m,n]t\k,l]Qi) — P(m; tj/c; 0)P(n; 0). This sum includes the 
contribution of paths which cross each other. In an annihilating random walk 
of otherwise non-interacting particles the contribution of all crossing paths 
have to be subtracted. Since we are on a infinite, one- dimensional lattice 
and both particles are identical this contribution is just the one given by all 
paths which start at site k and end at site n (instead of m) and which start 
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at site / and end at site m (instead of n). Therefore 

P(m, n; t\k, I; 0) = P(m; t\k; 0)P(n; t\l; 0) - P(n; t\k; 0)P(m; 0) (41) 

which is indeed what one obtains using the anticommutation relations in the 
free fermion approach. The same subtraction scheme generahzes to higher 
order conditional probabihties and is again conveniently captured in the free 
fermion anti-commutation relations. 

The contents of this section is the main result of the paper. In the follow- 
ing section we consider some of its consequences. Further relations can be 
obtained in a straightforward manner for multi-time correlation functions.^ 
The fact that H does not contain particle creation terms ensures that also 
the calculation m-point correlators for different times is reduced to the so- 
lution of the 2m-particle problem. The reduction of an equal-time m-point 
correlator to a 2m-particle problem for the one-dimensional model with ho- 
mogeneous hopping rates was first observed in using a mapping to a 



polymerization process. The new results presented here are the systematic 
and explicit way of expressing this reduction, its generalization to multi- 
time correlators and to non-translationally invariant systems and the final 
reduction to a single-particle problem in the dual environment. 

5 Density decay 

The result of the previous section is completely general as far as the hop- 
ping rates, observables and initial states are concerned. To be more spe- 

^In fact, all the results obtained here for equal-time correlations functions could have 
been phrased in more conventional language using the notion of duality as often done in 
the study of interacting particle systems. 
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cific we consider now as initial distribution an uncorrelated random initial 
state with density 1/2 in the sector of even particle number. This distri- 
bution is represented by the vector 1 1/2) = l/2^~-'^| sy"^"^ . We study the 
decay of the density in a so far unspecified, but fixed environment. One has 
( 1/2 \V = - i/2)^-^{0 I and thus 

p,{t) = ^{2\e'^'\k-l,k) (42) 

where ( 2 | = J2n>m ( ^ I is the sum over all states with two particles. 
Hence the density at site k at time t for a random initial distribution with 
density 1/2 is equal to one half of the survival probability of finding two 
particles anywhere in the dual environment, where initially two particles have 
been placed at sites k — l,k. The product form (^) of the two-particle 
conditional probability allows us to write 

Pkify 2 ^ ^ Pm,k ^Pm,k ~l~ Pm,k~l 

m 

oo 

^ ^ ^ ^ Pm,k Pm—n,k Pm—n—l,k—l ('^^) 
m 71=0 

with Pm,k = P{m] t\k] 0) as a short hand for the conditional probability. 

The expression (^Sf ) contains a double sum. This can be reduced to a 
single sum by considering the time derivative for the density for which one 

gets 

j^Pkit) = J](f„_i+C)(m-l,m|e-^*|A;-l,fc). (44) 

m 

The density at the boundary of a semi-infinite system, i.e. an infinite 
system with vanishing boundary hopping rates tq = £i = has an even 
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simpler expression, 

^ oo 

Pi{t) = -Y,Hn-A^,^)- (45) 

n=l 

In the dual system a particle which has moved to site cannot escape from 
there (except by annihilation with a second particle hopping from site 1 to 
site 0). Thus the boundary density at time t for a random initial distribution 
with density 1/2 is equal to one half of the survival probability of a single 
particle anywhere in the dual environment which has an absorbing boundary 
site and where initially the particle has been placed at sites 1. 

Eqs. (|42|) - (^5]) are exact and valid for any fixed environment. For further 
analysis we make two assumptions on the behaviour of the random walker. 
First consider site-symmetric processes with Vk = ik = Sk- In this case, the 
dual process is bond symmetric, ffc = ik+i, and therefore Pm,k = Pk,m- For 
environments such that limt_oo Zln Pm-n,k - Pm-n-i,k-i / Pm~n,k the 
second term in (^31) becomes small compared to the first for large times. This 
expresses the density at site k 

Pk{t) = i (P(A;; 2t\k; 0) + P{k; 2t\k - 1; 0)) . (46) 

in terms of return probabilities of a random walker. In the presence of an 
average drift with drift velocity v such that Pm,k = P2vt+k,m the analogue of 
(|6|) reads pk{t) = {P{k + 2vt; 2t\k] 0) + P{k + 2vt] 2t\k - 1; 0))/2. Here v 
may be itself be explicitly time-dependent. If furthermore to leading order 
in time P{k + 2vt; 2t\k] 0) = P{k + 2vt; 2t|/c — 1; 0), then one gets the simple 
result 

Pk{t) = P{k + 2vt;2t\k]0). (47) 
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The time-dependendence of the return probabihty depends on the environ- 
ment. 



6 Conclusions 

The main results of this paper are: 

(1) Diffusion-hmited pair annihilation and coagulation are equivalent to each 
other on any lattice and for arbritrary hopping and annihilation rates. The 
time-dependent density and density correlations of the mixed process have 
a simple expression in terms of the same quantities of the pure annihilation 
process (see Eq. (^BD). 

(2) In one dimension with nearest neighbour hopping and infinite annihila- 
tion rate the m-point correlation function for an arbitrary A^-particle initial 
state is given by correlation functions for a system with dual hopping rates 
(|36D with at most 2m particles initially. These correlation functions can be 
expressed in terms of single-particle conditional probabilities (see (P7|), (^Op), 



(41)). The derivation rests on a free fermion formulation of path integrals 
for annihilating random walks in one dimension. (3) Specifically for random 
initial conditions with density 1/2 one gets both the exact and approximate 
results of Sec. 5 the particle density. 

(4) For environments satisfying the assumptions detailed in Sec. 5 the den- 
sity at site k at time t with a random initial state is equal to the return 
probability of a single particle in the dual environment at time 2t (see (^^ ) . 

Expressing the density in terms of the return probability has been done for 
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computational convenience, not for physical reasons. Clearly, this relation is 
correct only in one dimension and only under the conditions outlined above. 
An open question which needs further investigation is the classification of 
hopping environments where these assumptions hold. It is also of interest 
to study universality of the density amplitude with respect to the initial 
density. This should be relatively straightforward in the approach developed 
above and is under investigation. A promising application of our results are 
disordered systems where further analysis requires taking a disorder 

average over products of single-particle conditional probabilities. Random 
walks in random environments are well-studied [^, 12| and, returning in 



our discussion to real systems such as exciton annihilation on polymers, one 
expects from this knowledge to gain considerable insight into the behaviour 
of diffusion-limited annihilation in one- dimensional disordered media. 

Acknowledgments 

The author would like to thank the Department of Physics, University of Ox- 
ford, where part of this work was done for kind hospitality and for stimulating 
discussions, particularly with J. Cardy and M.J.E. Richardson. 

References 

[1] Privman, V. (ed.): Nonequilibrium Statistical Mechanics in One Dimen- 
sion. Cambridge: Cambridge University Press, 1997 

[2] Lee, B.P.: J. Phys. A 27, 2633 (1994) 



25 



[3] 
[4] 

[5 



[7: 

[8 

[9 
[10 

[11 
[12 

[13: 



Schiitz, G.M., Mussawisade K., to be published 

Stinchcombe, R.B., Grynberg, M.D., Barma, M.: Phys. Rev. E 47, 4018 
(1993) 

Krebs, K., Pfannmiiller, M.P., Wehefritz, B., Hinrichsen, H.: J. Stat. 
Phys. 78, 1429 (1995) 

Schiitz, G.M.: J. Stat. Phys. 79, 243 (1995) 

Henkel, M., Orlandini, E., Schiitz, G.M.: J. Phys. A 28, 6335 (1995) 
Simon, H.: J. Phys. A 28, 6585 (1995) 
Peliti, L.: J. Phys. (France) 46, 1469 (1985) 
Cardy, J.: preprint |cond-mat 9607163] (1996) 



Alexander, S., Bernasconi, J., Schneider, W.R., Orbach, R.: Rev. Mod. 
Phys. 53, 175 (1981) 

Haus, J.W., Kehr, K.W.: Phys. Rep. 150, 265 (1987) 

Schiitz, G.M.: Integrable Reaction-Diffusion Processes and Quantum 
Spin Chains, to be published 



[14] Kadanoff, L.P., Swift, J.: Phys. Rev. 165, 310 (1968) 
[15] Doi, M.: J. Phys. A 9, 1465 (1976) 

[16] Grassberger, P., Scheunert, M.: Fortschr. Phys. 28, 547 (1980) 
[17] Racz, Z.: Phys. Rev. Lett. 55, 1707 (1985) 

26 



[18] Santos, J.E.: J. Phys. A, in press 

[19] Alcaraz, F.C., Droz, M., Henkel M., Rittenberg, V.: Ann. Phys. (New 
York) 230, 250 (1994) 

[20] Schiitz, G.M.: J. Phys. A 28, 3405 (1995) 

[21] Sponge, J.L.: Phys. Rev. Lett. 60, 871 (1988) 

[22] Cardy, J., Richardson, M.J.E.: in preparation 



27 



